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On  the  analysis  of  thermal  stresses  in  viscoelastic  solids* 

by 

Ell  Sternberg 
Brovm  University 


1 .  Introduction . 

The  study  of  thermal  stresses  and  defoimatlons  In  visco¬ 
elastic  solids  has,  for  various  technological  reasons,  attracted 
growing  attention  during  the  past  few  years .  This  paper,  which 
Is  devoted  exclusively  to  solids  that  under  isothermal  conditions 
and  for  Infinitesimal  strains  exhibit  linear  vlscoleastic  behavior, 
is  Intended  to  serve  a  dual  ptirpose.  First,  an  attempt  Is  made 
to  sketch  a  systematic  accowt  of  relevant  recent  theoretical 
developments;  second,  we  hope  to  supply  a  helpful.  If  necessarily 
Incomplete  guide  to  available  specific  results,  and  to  the 
literature  on  the  subject  at  hand. 

The  present  treatment  Is  largely  confined  to  the  quasi - 
static  analysis  of  thermal  stresses  in  homogeneous  suid  isotropic, 
linear  viscoelastic  media.  Questions  regarding  Inhomogeneous  and 
anisotropic  solids.  Inertia  effects,  and  thermo -mechanical  coupl¬ 
ing  effects,  are  touched  upon  only  at  the  end  of  the  paper 
(Section  7) . 


This  paper  was  prepared  xuider  Contract  Nonr-562(25)  of  Brown 
University  with  the  Office  of  Naval  Research  in  Washington,  D.C., 
for  the  International  Conference  on  High  Temperature  Structures 
and  Materials  at  Colvmbia  University,  New  York,  1963. 
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The  first  part  of  the  paper  (Sections  2,3t^)  deals  with 
thermo-vlscoelastlclty  theory  within  the  ass\imptlon  that  the 
mechanical  response  of  the  material  Is  temperature-independent. 
Since  the  rate  processes  of  viscoelasticity  are  known  to  be  highly 
sensitive  to  temperature  changes^  this  assumption  Is  remote  from 
physical  reality.  The  theory  of  temperature-independent  materials 
would  nevertheless  appear  to  merit  attention  as  a  useful  prelim¬ 
inary  to  a  more  realistic  treatment  of  the  problem  and  because  It 
has  In  fact  been  taken  as  the  basis  of  numerous  particular 
Investigations . 

Subsequent  portions  of  the  paper  (Sections  5>6)  deal 
with  t«iiperature-dependent  viscoelastic  solids,  particular  emphasis 
being  placed  on  the  theory  of  thermo-z*eologlcally  simple  materials. 
This  theory  rests  on  the  hypothesis  that  a  constant  temperature 
change  of  the  entire  solid  affects  Its  mechanical  response  merely 
within  a  \inlform  distortion  of  the  time  scale.  It  should  be 
apparent  that  the  smaller  amount  of  space  devoted  here  to  the 
Influence  of  temperature -dependent  response  characteristics  does 
not  reflect  the  relative  ln^)ortance  of  this  Influence  but  rather 
the  Increased  complexity  of  the  problem  once  the  temperature- 
dependence  of  the  material  Is  taken  Into  account. 

Although  throughout  the  following  discussions  some 
attention  Is  given  to  stress-strain  relations  In  differential- 
equation  form,  our  primary  preoccupation  la  with  the  more  general 
Integreil  laws  appropriate  to  solids  with  a  contlnuovis  spectrum  of 
relaxation  or  retardation  times.  Such  a  preference  Is  motivated 
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by  the  well-knovm  fact  that  the  finite  models  governed  by  the 
differential  operator  law  (despite  their  heuristic  V2ilue  and 
traditional  popularity)  provide  an  Inadequate  basis  for  the 
description  and  prediction  of  actual  viscoelastic  behavior  over 
any  appreciable  range  of  time  or  frequency.  This  shortcoming  Is 
even  more  pronounced  In  tanperature-dependent  models. 

The  following  treatment  of  thermo>vlscoelastlclty  theory 
Is  strongly  Influenced  by,  and  draws  heavily  on  the  results 
obtained  in,  an  earlier  study  of  the  Isothermal  theory  contained 
In  Reference  [1].  In  particular,  we  make  extensive  use  of  the 
properties  of  Stleltjes  convolutions  established  In  [l],  which 
enable  us  to  cope  economically  with  field  histories  that  are  dis¬ 
continuous  In  time.  As  emphasized  In  [l],  the  algebra  and 
calculus  of  convolutions furnish  the  natural  basic  tool  of  linear 
viscoelasticity  theory.  In  contrast,  an  excessive  reliance  on 
Integral  transforms  has  occasionally  tended  to  obsctire  the  Issues 
under  consideration.  Some  of  the  theoretical  results  Included  In 
this  paper,  thoiigh  elementary  in  character,  se«n  to  be  unavailable 
elsewhere.  In  this  sense  the  present  paper  is  not  wholly 
expository. 

2 .  Tonperature-lndependent  mechanical  response;  formulation  of 
boimdary-value  problqas. 

Hie  fundamented  system  of  field  equations  governing  the 
quasl-statlc  linear  theory  of  viscoelastic  solids,  in  the  absence 
of  thermo-mechanlcal  coupling  effects,  consists  of  the  linearized 
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displacement-strain  relations,  the  stress  equations  of  equilibrium, 

and  the  appropriate  linear  hereditary  stress -strain  law.  We  refer 

this  system  of  equations  to  rectangular  cartesian  coordinates 

2 

and  adopt  the  usual  Indlcial  notation. 

Let  u^(x,t),e^j(^t),(5^j(x,t)  be  the  component  values 
of  the  displacement.  Infinitesimal  strain,  and  stress  field 
histories  at  a  material  point  with  the  position  vector''  x  and  at 
the  time  t .  The  dlsplac«nent-straln  relations  then  become 


(2.1) 


whereas  the  stress  equations  of  equilibrium  assume  the  form 


0, 


(2.2) 


provided  denotes  the  ccmiponents  of  the  body-force  density  field 
history. 

With  a  view  toward  a  convenient  statonent  of  the  relevant 
constitutive  relations  for  Isotropic  viscoelastic  solids  we  Intro¬ 
duce  the  devlatorlc  components  of  strain  and  stress  by  means  of 


®1J  “  ®1J  ’  3  *lj®kk  '  ®13  *  ^IJ  "  3 
In  which  ^  Is  the  Kronecker  delta.  Next,  we  designate  by  T(x,t) 


* - 

Latin  subscripts,  \mless  otherwise  specified,  have  the  range  of  the 
Integers  (1,2,3);  summation  over  repeated  subscripts  Is  implied 
2md  subscripts  preceded  by  a  comma  Indicate  partial  differentia¬ 
tion  with  respect  to  the  corresponding  cartesian  coordinate. 

^Underlined  letters  designate  vectors. 
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the  local  Instantaneous  temperature  and  let  be  an  arbitrarily 
chosen  reference  temperature.  We  then  call  the  function  0 
defined  by 

©  =  T-T^  (2.4) 

the  temperature  field  history.  Finally,  we  write  o  for  the  coef¬ 
ficient  of  thermal  expansion,  which  we  shall  suppose  to  be 
constaint.  If  Instead  a  Is  a  fxmctlon  of  the  temperature,  all  of 
the  succeeding  considerations  undergo  an  entirely  elementary 
generalization:  In  this  Instance  one  need  merely  replace  o  by 
(2.4)  by 

fe>t) 

“  h  a(T«)dT'.  (2.5) 


The  appropriate  linear  Isotropic  hereditary  stress- 
strain  relations  In  the  form  of  the  relaxation  Integral  law  now 
appear  as 

“  ®1J*^^1  ^  ^klc  *  j^”3o0)*dG2 .  (2.6) 

Here  0^^  and  Qg  designate  the  respective  relaxation  moduli  In  shear 
and  Isotropic  compression.  These  moduli  are  at  present  fxmctlons 
of  the  time  exclusively  since  we  confine  our  attention  to  homo¬ 
geneous  solids  and  assiime  \mtll  later  on  that  the  material's 
mechanical  response  Is  temperature-independent.  In  writing  (2.6) 
we  have  employed  a  notation  for  Stleltjes  convolutions  Introduced 
previously  [l]  In  connection  with  the  Isothermal  theory.  Thus, 

If  f  and  g  are  functions  of  position  and  time,  u  =  f#dg  stands 
for  the  function  defined  by  the  Stleltjes  Integral 
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£‘(x,t-t‘)dg(xit«),  (2.7) 


t  *=-00 

provided  this  Integral  is  meanlngf tjl . 

We  observe  that  according  to  (2.6),  the  local  stress 
tensor  at  each  fixed  Instant  is  a  linear,  continuous,  and  isotropic 
functional  of  the  entire  preceding  local  strain  and  temperature 
histories;  moreover,  this  mapping  of  the  strain  and  temperature 
hsltorles  Into  an  associated  history  of  stress  is  invariant  under 
a  translation  of  the  time  scale  and  has  the  property  that  the 
stresses  Induced  by  a  free  thermal  expansion  vanish  identically. 
One  can  show^  (tinder  suitable  regularity  assumptions)  that  (2.6) 
is  the  most  general  stress-strain  law  conforming  to  the  foregoing 
requirements . 

The  field  equations  (2.1), (2. 2), (2. 6)  must  hold  through¬ 
out  the  space-time  domain^  Rx (- oo,  co),  l.e.  for  all  (x,t)  such 
that  X  is  in  the  open  region  R  occupied  by  the  interior  of  the 
body  and  t  lies  in  the  interval  (-00,00).  To  these  field  equa¬ 
tions  we  may,  without  essential  loss  of  generality,  adjoin  the 
requirement  that  the  body  is  originally  imdisturbed  in  the  sense 
of  the  initial  conditions 


u-  =  0  =  0  on  Rx(-oo,o),  (2.8) 

1  ij  IJ  1 


See  [1],  Section  2,  for  the  corresponding  result  in  the  iso- 
theiroal  theory. 

^  We  use  the  conventional  notation  A  x  B  for  the  cartesian  product 
of  a  set  A  and  a  set  B. 
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where  R  stands  for  the  closure  of  Ci,  l.e.  for  the  union  of  R  and 
Its  boundary.  Finally,  in  the  case  of  the  standard  mixed  bovindary- 
value  problem,  the  bovmdary  conditions  become 

Ui  »  u^  on  B^x(- 00,00), 

Si  »  on  B2x(-oo,oo). 

Here  and  Bg  are  complementary  subsets  of  the  boundary  B  of  18 , 

■w 

nj  is  the  outward  xanit  normal  of  B,  whereas  u°  and  are  pre¬ 
scribed  surface  displacements  and  surface  tractions,  respectively. 

In  the  absence  of  an  explicit  statement  to  the  contrary, 
the  point  sets  R,B^,B2  will  henceforth  be  asemned  to  be  independ¬ 
ent  of  the  time.  Further,  we  shall  suppose  from  here  on  that  R 
U  a  .ounae.  ««aXar  o.  apace®  ana  tna*  ane  Potn 

Integrable . 

The  problem  \inder  consideration  thus  consists  in  deter¬ 
mining  field  histories  which  —  for  given  R,B^,Bg, 

known  Oj^,Og,a,  and  prescribed  F^,0,u^,S^  —  satisfy  the  field 
equations  (2.l),(2.2),(2.3)a(2.6)  in  {2x(-oo,  00)  and  meet  the 
initial  conditions  (2.8),  as  well  as  the  boundary  conditions  (2.9). 
The  t«nperature  field  history  0  may,  in  partlcvilar,  be  specified 
as  the  solution  of  an  Independent  heat-conduction  problem. 

For  future  economy  we  now  adopt  the  following  definition 

of  a 

^y  a  "regular  region  of  space"  we  mean  a  region  whose  boundary 
consists  of  a  finite  ntimber  of  non-intersecting  closed  regvilar 
surfaces",  the  latter  term  being  used  in  the  sense  of  Kellogg  [2]. 
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Viscoelastic  state .  We  say  that  the  ordered  array  of  field 
histories  J  belongs  to  the  class  of  viscoelastic  states 

on  Rx(-oo,  oo)  corresponding  to  the  data  and  write 

[Ui,eij,dij]€YtV^2'“'^l'®^  on  Rx(_oo,  oo)  (2.10) 

(a)  Qp(p=l,2)  vanishes  on  (-00,0),  Is  twice  continuously 
differentiable  on  [O,oo),  and  Qp(0)  >  0; 

(b)  vanish  on  Rx(-oo,o)  and  are  continuous  on 
Rx[o,oo)^  0  being  once  —  and  u^^  three  times  continuously 
differentiable  on  Rx[0, oo)* 

(c)  equations  (2.l),(2.2),(2,3)*(2,6)  hold  In  Qx(-oo,oo). 

If.  in  particular.  0*0  on  Rx(-oo^  °°)>  we  say  that  the  state 
[Ui,eij,<Jij]  Is  Isothermal  and  write 

fuj^^e j  3  ^ ^*^1  ^  Rx( -00^  00) .  (2 .11 ) 

Requirements  (a),{b),(c)  are  evidently  partly  red\andant 
but  mutually  consistent.  Also^  the  smoothness  assumptions  con¬ 
tained  In  (a),(b)  could  be  relaxed  (especially  as  far  as  the  time- 
dependence  Is  concerned)  at  the  expense  of  more  elaborate  regular¬ 
ity  hypotheses.  Such  refinements  would,  moreover,  necessitate 
occasional  supplementary  smoothness  assmnptlons  In  subsequent 
theorans  and  would  merely  tend  to  detract  from  the  main  purpose 
of  the  present  paper. 

It  Is  essential  to  observe  that  the  field  histories 
Involved  In  the  preceding  definition  of  a  viscoelastic  state  may 
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evidently  exhibit  finite  Jmp  discontinuities  at  t=0.  Although 
such  singular  field  histories  are  physically  unrealistic  — 
partlculsu?ly  within  the  context  of  the  quasi-static  theory  —  they 
play  an  Important  part  In  the  theory  of  Integration  of  the  govern¬ 
ing  field  equations  because  of  Duhamel’s  principle.  One  of  the 
advant8iges  derived  from  the  formulation  of  the  constitutive 
relations  (2.6)  In  terms  of  Stleltjes  convolutions  lies  In  the 
fact  that  It  permits  a  systematic  treatment  of  the  relevant  dls- 
contlniiltles  and  enables  one  to  avoid  the  usual  purely  formal 
manipulations  with  Dirac *s  delta  function. 

If  (2.10)  holds,  the  stress-strain  law  (2.6)  admits  the 
conventional  Rl«nann  Integral  representation  (cf.  [l],  Oheorem  3.^)# 
valid  for  (x,t)  In  Rx[0,oo), 

8ij(x,t)  «  03^(t)e^j(x)  +  I  03L(t-t')ej^j(x,t»)dt», 

0 

<Jjjjj(x,t)  »  ’  3a0(x)] 

'  f 

+  J  Q2(t-t*)[ejjjj.(x,t«)  -  3o®(x,t‘)]dt». 

0 

Here,  as  In  the  sequel,  f  stands  for  the  first  time  derivative  of 
a  function  f  of  position  and  time,  whereas 

f(x)  “  f(x,0) .  (2.13) 

We  consider  next  the  case  of  an  elastic  solid.  To  this 
end  let  h  henceforth  denote  the  Heaviside  \anlt  step  functi.on 
defined  by 
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h(t)  =  0  for  -  oo  <  t  <  0, 
h(t)  =  1  for  0  <  t  <  oo. 

It  Is  Immediate  from  (2.12)  that  for  the  particular  choice 


(2.14) 


0^  =  2nh  ,  Og  =  3J<h  ,  (2.15) 

where  (j.  and  k  are  constants^  the  stress-strain  relations  (2.6) 
pass  over  Into 

®lj  ~  *  ^kk  ~  3Jt(Ej^”3o0)  •  (2.16) 


Equations  (2.l6)j  In  view  of  (2.3)>  are  equivalent  to  Hooke’s  law 
(modified  to  accoxint  for  thermal  expansion)  j  provided  p.  and  x  are 
the  shear  modtilus  and  the  hulk  modulus  of  the  elastic  material^ 
respectively.  This  remark  reveals  the  role  of  the  classical 
quasi -static  theory  of  thermoelasticity  as  a  special  case  of  the 
theory  of  thermo-vlscoelasticlty  laider  consideration  and  motivates 
the  following  definition. 

Quasi -static  elastic  state .  If  (2.10)  holds  emd  the  relaxation 
functions  obey  (2 .15),  In  which  u  and  x  are  (positive)  constants . 
we  say  that  belongs  to  the  class  of  quasl-statlc 

elastic  states  on  Rx(-oo, 00)  corresponding  to  the  data  p,x,a,P^,0, 
and  write 


[Ui,eij,qij]^g  [|i,H,a,Pi,0]  ~  Rx(-oo,oo) 
Q 


(2.17) 


It  Is  clear  fron  the  above  definition  that  a  quasl- 
statlc  elastic  state  which  falls  to  vanish  Identlcedly  on 
Rx(-oo^  00),  cannot  be  Independent  of  time  on  this  entire  space- 
time  domain.  For,  (2.17)  Inqjlles  that  ®wst  meet  the 
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Initial  conditions  (2.8).  Since  we  shall  need  to  deal  also  with 
purely  position-dependent  solutions  of  the  fundamental  field  equa¬ 
tions  of  thennoelastlclty  theory,  we  find  It  convenient  to  Intro¬ 
duce  the  notion  of  a 

Stationary  elastic  state.  We  say  that  ]  belongs  to 

the  class  of  stationary  elastic  states  on  R  corresponding  to  the 
da^  |i,x,a,P^,0,  and  write 

^^1'®1J'^1J^^  ^  ^  (2.18) 

if  J 

(a)  p.  and  x  are  (not  necessarily  real)  constants; 

(b)  are  functions  of  position  continuous  on 
R,  0  being  once  —  and  u^  three  times  continuously  differ¬ 
entiable  on  R; 

(c)  equations  (2.l),(2,2),(2,3)#(2.l6)  hold  In  fS  . 

Ohe  reason  why  we  do  not  restrict  and  x  to  be  real 
constants  In  the  present  Instance  will  become  apparent  at  the  end 
of  Section  4.  Meanwhile  we  note  that  (2.1?)  Implies,  for  each 
fixed  t  In  (-00,00), 

[UjL(-,t),e^j(*,t),<J^j(*,t)]€£[^,H,a,PjL(‘,t),e(*,t)]  on  R  (2.19)^ 

so  that  every  quasl-statlc  elastic  state  on  Rx(-oo,  00)  may  be 
regarded  as  a  one-parameter  family  (with  time  as  the  parameter) 
of  stationary  elastic  states  on  R  corresponding  to  the  same 

^  If  f  is  a  function  of  position  and  time  defined  on  Rx(-oo,  00), 
we  write  f(*#t)  for  the  fxmctlon  of  position  defined  on  R  that 
resvilts  from  the  mapping  f  by  holding  the  time  fixed  in  (-00,00) 
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elastic  constants  and  to  the  appropriate  family  of  body-force  and 
temperature  distributions.  This  observation  reflects  the  fact  that 
the  time  In  quasl-statlc  elasticity  theory  plays  the  role  of  a  mere 
parameter,  whereas  the  quasl-statlc  theory  of  viscoelasticity  Is 
genuinely  time-dependent. 

We  have  so  far  based  the  fomiulatlon  of  boundary-value 
problems  In  thermo-vlscoelasticlty  theory  upon  the  relaxation 
Integral  law  (2.6).  On  the  assumption  that  Qp(p=l,2),  and 

meet  the  conditions  set  down  under  (a),(b)  In  the  definition 
of  a  viscoelastic  state,  the  stress-strain  relations  (2.6)  may 
be  inverted.  This  leads  ([1],  Theorem  3.3)  to  the  equivalent 
creep  Integral  law 

*  (2 .20) 

where  and  Jg  are  the  creep  compliances  in  shear  and  isotropic 
compression,  respectively.  Moreover,  the  two  pairs  of  response 
functions  0^  and  Jp  (P=l,2)  are  linked  by  the  relation 

Gp*dJp  =  h  on  (-00,00)  (2.21) 

or,  using  a  notation  adopted  In  [1]  (Hieorem  1.3)  for  the 
"Stleltjes  Inverse"  of  a  function  of  time, 

Jp  =  Gp^  on  (-00,00).  (2.22) 

Eqiiatlon  (2 .21)  Implies  further  that 

f  • 

GpJp(t)  +  J  Gp(t-t')Jp(t»)dt«  =1  for  0  <  t  <  00.  (2.23) 

0 

Bie  familiar  physical  significance  of  the  relaxation 
moduli  and  of  the  creep  compliances  Is  Immediate  from  (2.12)  and 
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Its  creep  counterpart.  Thus 

®1J  “  ^  Rx(-oo,oo)  Implies  s^j  *  Gr^  on  Rx( -00,00),  (2.24) 
®1J  “  ^  Rx(-oo,oo)  Implies  e^^j  “  on  Rx(-oo,  00).  (2.25) 

Analogous  Interpretations  apply  to  Qg  and  Jg.  If  f  Is  a  function 
of  position  and  time  that  possesses  a  Laplace  transform  with 
respect  to  time,  we  write 

00 

f(xjT))  =  Jt{f(x,t) ;tj}  *  I  f (2^t)exp(-Tit)dt,  (2.26) 


in  which  T|  designates  the  transform  parameter.  Accordingly, 
supposing  G|p  and  Jp  (P«®1,2)  to  he  of  exponential  order  as  t  —^00^ 
one  draws  from  (2.23)  in  conjtinctlon  with  the  convolution  theorem 
for  the  Laplace  transform  that 


^p(4)7p(Ti)  = .  (M,2).  (2.27) 

We  tvim  finally  to  stress-strain  relations  In  differ¬ 
ential  equation  form,  i.e.  to  what  is  traditionally  referred  to 
as  the  differential  operator  law  of  linear  viscoelasticity.  For 
this  purpose  we  first  make  the  following  notational  agreements. 

If  f  Is  a  (suitably  smooth)  function  of  position  and  time  defined 
on  Rx(-oo,oo),  we  denote  Its  n-th  partial  time  derivative  by 

f(n)  ^  jjnf  (n=0,l,2, ...),  (2.28) 

0/  \ 

D  being  the  time-derivative  operator.  Further,  we  write  f'  '  for 
the  function  of  position  defined  by 


f^^'^Cx)  =  f^*'^(x,0+-)  for  X  In  R. 


(2.29) 
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Hie  differential  operator  law  now  becomes 


Pi(D)8^j  "  •  r2(»)»lclc  ■  «2(»)[8ac-3ael,  (2.30) 

provided  Pp(D),Qp(D)  (p=l,2)  are  the  linear  differential  operators 

Np  Np 

“  n=o  “  n-o 

Here  Np  (p*l,2)  is  a  non-negative  integer  and  the  coefficients 
pp.n^Qp.n  constant  response  parameters  for  a  given  (tempera- 
tiire- Independent)  material.  Moreover,  we  may  evidently  assume 
that  either  P^.p  ^  0  ^.p  0  when  n  «=  Np  (p=l,2)  so  that,  for 

fixed  p,  at  least  one  of  the  operators  in  (2. 31)  has  the  degree 
Np  .  The  differential  eq\2ation8  (2.30)  are  to  be  met  in  Rx(0,  oo) 
and  must  be  accompanied  by  the  following  initial  conditions. 


which  are  valid  on  R: 


n-r 

^2  0/  \ 
Z  p„ 

in  kk 


2  q, . 
n*»r 


*l;n''lj 


(1*1,2,  . .  .,N^), 


N 

'  3a0('^‘"’^Kr=l,2,...,N2). 

ns*r  * 


>(2.32) 


As  is  clear  from  an  elementary  extension  of  Theorem  4.1 
in  [1],  equations  (2 .30), (2 .32)  are  Implied  by  (2.8)  and  the 
relaxation  integral  law  (2.6)  or  the  creep  Integral  law  (2.20) 
provided  sufficiently  smooth  £uid  according  as 

the  relaxation  moduli  or  the  creep  compliances  exhibit  the  de¬ 
generacies  characteristic  of  a  finite  spectrum  of  relaxation  or 
retardation  times.  In  the  first  case  p^.p  /  0  for  n*Np  (P»l,2), 
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whereas  In  the  second  Instance  /  0  for  n  -  Np  (p«l,2).  The 
physical  significance  of  the  Initial  conditions  (2.32)  was 

g 

established  In  [l](Theorem  4.2)  for  the  isothermal  case  .  Thus 
suppose  @  veinlshes  Identically  and  hold  x  fixed.  Then^  rovighly 
speaking,  (2.32)  are  necessary  and  sufficient  that  every  pair  of 
local  histories  e^j(x,  •)#  which  vanishes  on  (-oo,o), 

satisfies  (2.30)  on  (O,  00),  and  exhibits  finite  Jianp  discontinu¬ 
ities  at  t«0,  be  the  limit  of  a  pair  of  history  sequences  which 
vanishes  on  (-00,0),  obeys  the  same  stress-strain  relations  (2.30), 
and  has  time  derivatives  of  the  orders  entering  (2.30)  that  are 
continuous  on  ( -  00,  00) . 

If  an  Integral  law  (2.6)  is  reducible  to  a  differential 
operator  law  (2.30),  then  the  Laplace  transforms  of  the  relaxation 
fvinctions  exist,  are  necessarily  rational,  and  are  given  by  ([1], 
Theorem  4.8) 

Qp(n) 


(2.33) 


Analogously,  the  reduclblllty  of  a  creep  IntegrsJ.  law  (2 .20) 
implies 


p 


(2.34) 


Additional  results  concerning  the  transition  from  2in  Integral  to 
a  differential  operator  law,  and  vice  versa,  may  be  found  in  [l] 
(Section  4) . 


F 


See  also  Boley  and  Weiner  [3],  Art.  15*6,  for  a  related  discus¬ 
sion. 
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Since  materials  with  a  finite  releucatlon  or  retardation 
spectrum  admit  the  familiar  model  representation  In  terms  of 
easily  visualized  finite  networks  of  springs  and  dashpots,  the 
bulk  of  the  existing  literature  on  viscoelasticity  has  favored 
the  differential  operator  law  over  stress -strain  relations  In 
Integral,  form.  It  Is  well  to  bear  In  mind,  however,  that  such 
a  preference  Is  not  warranted  on  the  grounds  of  theoretical 
generality;  nor  Is  It  Justified  from  a  practical  point  of  view 
since  an  adequate  description  of  the  behavior  of  actual  visco¬ 
elastic  solids  over  any  appreciable  range  of  time  (or  frequency) 
ordinarily  necessitates  the  use  of  differential  operators  of  a 
comparatively  high  order. 

3.  Temperature -Independent  mechanical  response;  general  Implica¬ 
tions  . 

We  turn  next  to  general  results  concerning  the  funda¬ 
mental  field  equations  and  boundary-value  problems  discussed  In 
the  preceding  section.  In  this  connection  we  shall  limit  our 
attention  to  materials  governed  by  the  relaxation  Integral  law 
(2.6);  anadogous  conclusions  apply  to  the  creep  Integral  law 
(2.20)  and  to  the  differential  operator  law  (2.30).  To  facilitate 
the  task  at  hand  we  first  establish  a  link  between  the  theory  of 
thermo-vlscoelastlclty  under  present  consideration  and  the  corre¬ 
sponding  Isotheiroal  theory,  which  Is  supplied  by  the 
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Body-force  analogy .  Suppose 

on  Rx(-oo,  oo)  (3.1) 

and 

=  u“  on  B^x(-oo,  oo),  =  S°  ^  B2x(-  oo,  oo) .  (3.2) 


Further j  let  ^  field  histories  defined  on  Rx(-oo,oo 

by  means  of 


Uf  *  Ui  ,  =  e^j  , 


"ij  '  '  ^1  ‘  ^1  -  “®,l«d02 


y  (3.3) 


Then 


and 


[uiieij,3ij]6YCV®2'^l^  on  Rx(-oo,oo)  (3.4) 


Uj^  =  u^  on  B^x(-oo,  oo),  +  on^^S^dOg  on  BgX(-oo,  oo).  (3.5) 


Conversely,  If  0  has  the  same  properties  as  In  the 
definition  of  a  viscoelastic  state,  then  (3.3)#(3.4),(3.5)  Imply 

(3.1)  (3.2). 

The  truth  of  the  foregoing  theorem  Is  readily  Inferred 
from  the  underlying  definitions  of  viscoelastic  and  Isothermal 
viscoelastic  states  with  the  aid  of  a  known  result  ( [1 ], 

Theoran  1.6)  on  the  space  differentiation  of  Stleltjes  convolu¬ 
tions.  The  theorem  yields  as  a  special  case  the  well-known  body- 
force  analogy  of  thermoelastlclty^  If  the  relaxation  functions 
obey  (2.15).  Although  the  generalized  analogy  xmder  discussion. 


9 


See,  for  example,  [3]j  Art.  3.3. 
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which  reduces  the  mixed  boundary-value  problem  of  thermo-visco¬ 
elasticity  to  the  standard  mixed  problem  in  the  isothermal  theory, 
is  of  no  practical  usef\ilness  as  far  as  the  actual  solution  of 
specific  problems  is  concerned,  its  theoretical  Interest  is  con¬ 
siderable.  For,  the  analogy  enables  one  to  obtain  economically 
extensions  to  non-lsothermal  conditions  of  theorems  available  in 
the  Isothermal  theory.  We  now  cite  some  of  the  more  Important 
conclusions  reached  in  this  manner. 

Volterra’s  [4]  theorem^^  concemlrig  the  uniqueness  of 
the  solution  to  the  mixed  Isotheiroal  problem,  in  conjunction  with 
the  body-force  analogy,  at  once  furnishes  the 
Uniqueness  theorem .  Suppose 


^ j > (3 ^ j 3 Qg, c, 3  8x(  —  oo,  oo), 

~  Rx(-oo,  oo), 

and  let 

Uj^  “  u^  on  B^x(-oo,  oo),  on  B2x(-oo,  oo) 


r(3.6) 


(3.7) 


Then 


[ui,eij,<3ij] 


^^I'^lJ'^^lj^  Rx(-oo,  oo),  (3.8) 


where  w^  =  0  Rx(-oo,o)  and  w^^  represents  an  ( infinitesimal ) 
rigid  motion  of  the  entire  body  on  Rx[0, oo) . 

According  to  (a)  in  the  definition  of  a  viscoelastic 
state,  the  relaxation  functions  entering  (3.6)  are  reqxalred  to 


See  [l].  Section  8,  for  a  slightly  more  general  version  and  a 
more  detailed  proof  of  Volterra's  result  in  the  special  case 
of  an  isotropic  solid. 
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satisfy  the  inequalities 

Op  >  0  (P»l,2),  (3.9) 

l.e.  must  possess  positive  Initial  values.  The  striking  fact 
that  no  restrictions  (apart  from  a  smoothness  assumption)  need 
be  Imposed  on  the  subsequent  behavior  of  Op  in  order  to  lns\ire 
uniqueness^  will  be  made  plausible  shortly. 

A  characterization  of  the  initial  field  distributions 
appropriate  to  a  viscoelastic  state  is  supplied  by  the  following 
theorem,  which  is  an  elementary  consequence  of  Theorem  6.1  in  [l] 
and  of  the  body-force  analogy. 

Initial  response.  Let  3  meet  (2.10).  Then 

0  0  0  o  o 

[Ui,eij,d^j3  e&[li,x,a,Pi,0]  on  R  (3.10) 


with 


(3.11) 


This  result  lends  a  precise  meaning  —  within  the 
context  of  thermo-viscoelasticity  theory  —  to  the  familiar  asser¬ 
tion  that  the  initial  response  of  a  viscoelastic  solid  is  elastic . 
Moreover,  the  preceding  theorem  enables  one  to  determine  directly 
the  initial  displacements,  stresses,  and  strains  belonging  to  the 
solution  of  a  mixed  boundary-value  problem  in  thermo-viscoelastic - 
Ity  theory  by  solving  a  steady-state  theraoelastlc  problem;  the 
latter  is  governed  by  the  initial  body-force,  temperatrire,  and 
stirface  data  of  the  original  problem,  as  well  as  by  the  elastic 
constants  (3.11).  Similarly,  Theorem  6.2  in  [1]  furnishes  an 
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analogous  characterization  of  the  Initial  (right-hand)  time 
derivatives  of  all  existing  orders  belonging  to  the  desired  visco¬ 
elastic  state.  These  Initial  time  rates  (n«l,2,...) 

may  be  found  directly  from  the  corresponding  Initial  data  by 
solving  successively  a  sequence  of  steady-state  thermoelastic 
problems,  each  of  which  Is  governed  once  again  by  the  elastic 
constants  defined  In  (3.11).  Since,  according  to  the  uniqueness 
theorem  of  the  classical  (uncoupled)  theory  of  thermoelasticity, 
the  Inequalities  p.  >  0,  x  >  0  ai*e  sufficient  to  guarantee  the 
vinlqueness  of  a  suitably  regxilar  solution  to  the  mixed  problem, 
Volterra‘s  uniqueness  criterion  (3.9)  Is  not  so  surprising. 

An  additional  conclusion  regarding  the  time-dependence 
of  viscoelastic  states,  that  permits  one  to  Infer  the  smoothness 
of  the  state  with  respect  to  time  from  the  corresponding  smooth¬ 
ness  of  the  field  and  surface  data.  Is  deduclble  from  Theorem  6.4 
In  [l].  We  proceed  now  to  a  theorem  on  the 
Position-dependence  of  viscoelastic  states.  Let 
meet  (2.10)  and  assvime  that 


V.p  =  0,  VaP  «  0,  7^0  *  0  on  i2j<-oo,  oo)  . 

(3.12)^^ 

Then  throughout  fi)x(-ooioo). 

7^(7 -u)  *  0,  7^(7-u)  =  0, 

(3.13) 

’V  "  ’'"13  *  ’'”13  ■ 

(3.1^) 

TT - 

Here,  as  In  the  sequel,  V  Is  the  usual  spatial  gradient 

operator,  whereas  7*,  V*,  and  ^  designate  the  divergence, 
curl,  and  Laplaclan  operators,  respectively. 
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Ihus  in  the  present  circvmiatances,  as  In  thermoelastielty 
theory,  (3.12)  Imply  that  the  dilatation  and  the  rotation  are 
harmonic,  while  the  cartesian  components  of  displacement,  strain, 
and  stress  are  blharmonlc  functions  of  position.  It  follows  that 
all  of  these  field  histories  possess  continuous  space  derivatives 
of  all  orders.  Ihe  above  theorem  is  a  consequence  of  Theorem  6.6 
in  [1]. 


The  comterpart  of  Betti  ‘  s  elastostatlc  reciprocal  rela 
tions  in  the  isothermal  theory  of  viscoelasticity  given  in  [1] 
(Theorem  7.^)*  in  view  of  tbe body-force  analogy,  the  divergence 
theorem,  and  by  virtue  of  Theorems  1.2,  1.6  in  [l],  leads  to  the 
following 

Reciprocal  theorem .  Suppose 


[Ui,eij,dij]  €Y[a^,CI^,a,P^,0]  on  Rx(-oo,oo), 
[Uiie^j,d^j]€T[C^,O2,o,P^,0‘ ]  on  Rx(-oo,oo). 


r (3.15) 


Then,  on 


t 


(-00,  oo), 

p 

S,*duJdA  +  P,»du,*dV  +  o  , 

^  R  ^  R 

J  S%du.dA  +  P*#du,dV  +  o  , 
B  ^  R  ^  R 

( 


@«de 


0»deii»dOgdV  = 
0**de^j^*dQ2dV  = 
»da2dV  « 


J 


0  #de^^*dQ2dV. 


(3.16) 
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It  should  be  emphasized  that  the  reciprocal  relations 

obtained  from  the  present  theorem  by  specializing  Op  (P=l,2)  In 

accordance  with  (2.15)  and  by  setting  0“O  on  Rx(- 00,00),  are 

distinct  from  the  quasl-statlc  adaptation  of  Betti's  reciprocal 

relations  In  Isothermal  elastostatlcs .  This  distinction  disappears 

If  the  body  force  and  temperature  field  histories,  as  well  as  the 

surface  data,  of  the  two  states  In  (3*15)  are  separable  functions 

of  position  and  time  that  share  a  common  time -dependence.  In 

this  particular  Instance  the  convolutions  entering  (3.16)  (e.g. 

1 .  12 

S^*dUj|^J  can  be  shown  to  be  reducible  to  ordinary  products 

(e.g.  S^u^),  and  the  resulting  reciprocal  relations  are  Identical 

with  those  arising  from  the  extension  of  Betti's  theorem  to  quasl- 

statlc  thermoelasticity  theory. 

The  reciprocal  theorem  stated  above  may  be  used  to 

derive  formulas  for  the  average  strains  Induced  In  a  (temperature- 

independent)  viscoelastic  solid  by  given  body-force  amd  tempera- 

10 

ture  field  histories  and  by  specified  surface  tractions.-'  We 
Include  here  merely  a  formula  for  the  total  volume  change,  which 
can  be  deduced  directly  from  Uieorem  7.6  In  [l]  with  the  aid  of 
the  body-force  analogy  and  the  divergence  theorem. 

Cf.  theorem  7-5  In  [l]. 

See  [5]  for  a  derivation  of  suialogous  resvilts  In  the  Isothermsd 
theory  of  anisotropic  viscoelastic  solids. 
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Volume  change.  Prom  (2.10)  follows ^  for  t  ^  (-00,  00), 

AV(t)  «  [A*dJ2](t)  +  3a  I  ®(x,t)dV,  (3.1?) 

provided  Jg  the  creep  compliance  corresponding  to  the  relaxa¬ 
tion  modulus  Qg  and  where 

AV(t)  =  J  e..(x,t)dV,  (3.18) 

R 

A(t)  «  I  x^Sj^(x,t)dA  +  I  Xj^P^(x,t)dV.  (3.19) 

It  Is  essential  to  observe  that  (3.1?)  niay  also  be 
derlved^^  solely  from  (3.l8),(2.2),  and  the  second  of  (2.6)  by 
recourse  to  (2 .21)  and  Theorems  1.2,  1.6  In  [1].  Consequently, 
formula  (3.17)  Is.  valid  for  any  solid  that  is  In  a  state  of 
Infinitesimal  deformations  and  has  a  temperature-independent 
linearly  viscoelastic  dllatatlonal  response,  regardless  of  the 
nature  of  Its  response  In  shear.  If,  In  particular,  Sj^=0  on 
3x(-oo,oo)  and  P^=0  on  Rx(- 00,00),  (3.17)  Implies 

AV(t)  “  3a  J  0(x,t)dV  for  t  In  (-00,00).  (3-20) 

R 

The  same  conclusion  was  also  reached  by  Nowackl  [7](Chapter  XI)  on 
the  basis  of  a  less  direct  argument.  Thus,  In  the  absence  of 
surface  tractions  and  body  forces,  the  total  vol\ame  change  Is  the 
same  as  that  which  would  arise  from  uninhibited  thermal  expansion. 

See  [6 3,  where  such  an  alternative  derivation  Is  given  for  the 
special  case  of  a  material  with  a  purely  elastic  dllatatlonal 
response . 
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15 

This  well-knovm  result  of  classical  thermoelastlolty  theory  ^  Is 
accordingly  seen  to  survive  without  modification  under  the  present 
more  general  hypotheses . 

The  next  theorem^  which  also  pertains  to  vanishing 
surface  tractions  and  body  forces^  follows  from  the  appropriate 
specialization  of  a  result  established  In  [6]  for  a  class  of 
temperature -dependent  viscoelastic  materials. 

Stress -free  temperature  field  histories .  Assume  (2 .10)  holds  and 
suppose  on  Bx[o,  oo),  R^[0,oo) .  Then 

Rx(  -  oo,  oo)  ^  and  only  If 

0(x,t)  =  +  aj^(t)xj^  for  (x,t)  3n  Rx[0,oo),  (3.2l) 

where  a^^a^  are  functions  continuous  on  [0, oo) . 

An  Indication  of  the  proof  of  the  necessity  of  (3.2l) 

If  the  body  Is  to  be  free  from  stress,  appeared  prior  to  [6]  In 
a  note  by  Hilton  [9].  Analogous  two-dimensional  theorems  on 
temperature  distributions  that  fall  to  Induce  thermal  stresses 
(applicable  to  conditions  of  plane  strain  and  generalized  plane 
stress)  are  also  to  be  found  In  [6].  For  the  familiar  correspond¬ 
ing  three-dimensional  and  two-dimensional  results  In  thermo - 
elasticity  theory,  which  are  Implied  by  the  theorems  under  present 
consideration,  reference  may  be  made  to  [33#  Arts.  3.9#  4.9. 

Before  concluding  this  section  we  merely  mention  several 
additional  Implications  of  the  body -force  analogy.  Thus,  one  may 

See,  for  example,  [3]#  Art.  9.15.  The  result  was  reached  by 
Hleke  [8],  who  based  his  proof  on  an  extension  to  thermo- 
elastlclty  theory  of  Betti’s  reciprocal  relations. 
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extend  Duhamel’s  principle  of  the  isothermal  theory  (see  [l]. 
Section  3)  to  the  analysis  of  thermal  stresses  and  deformations. 
Similarly^  one  may  obtain  analogous  extensions  of  the  variational 
principles  of  isothermal  viscoelasticity  theory Finally,  the 
Integral  representations  (in  terms  of  Green's  functions)  of  the 
solution  to  Isothermal  boundary-value  problems  established  in 
[11]  are  readily  generalized  to  include  temperature  effects  within 
the  current  theoretical  framework. 

Some  further  consequences  of  the  body-force  analogy 
will  be  encountered  in  the  succeeding  section,  which  deals  with 
results  that  have  a  direct  bearing  on  the  solution  of  boundary- 
value  problems  of  the  type  formulated  in  Section  2. 

4.  Temperature-independent  mechanlcetl  response;  methods  of 
integration «  available  solutions. 

With  a  view  toward  alternative  formulations  of  the 
bo\andary-value  problems  stated  in  Section  2,  we  now  cite  the 
appropriate  displacement  equations  of  equilibrium  and  stress 
equations  of  compatibility,  which  could  be  arrived  at  by  a 
respective  elimination  of  either  the  stresses  and  strains,  or  of 
the  displacements  and  strains,  among  the  governing  fundamental 
system  of  field  equations.  Actually,  these  desired  generalizations 
of  the  Cauchy  equations  of  equilibrium  and  of  the  Beltrami -Mlchell 
eqviatlons  of  compatibility  are  immediate  from  their  isothermal 
counterparts  in  Theorem  5.5  and  Theorem  5.7  of  [l]  because  of  the 
body-force  analogy. 

See  Ourtln  [10]  for  a  generalization  to  Isothermal  viscoelastic¬ 
ity  theory  of  the  classical,  and  of  more  recent,  alastostatlc 
variational  principles. 
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Displacement  equations  of  equilibrium.  Prom  (2.10)  follows j  on 
®x(-  oo,  oo), 

or,  ea\ilvalently« 


V^u*dQj^  +  ^  VV*u*d( 0^+202)  +  2P  «  2o70iKaa2-  (^*2) 

Stress  equations  of  compatibility .  Prom  (2.10)  follows j  on 
fJx(-oo,  00), 


(4.3) 


where  Jp  Is  the  creep  compliance  corresponding  to  Op  (p*l,2), 

-  a»^.7®eiKia2»a(!Hjj^)  -  ae^ij» 

and 

Q  »  .  (4. 5)^"^ 


The  mixed  problem^  when  cast  In  terms  of  displacements 
alone,  consists  In  finding  a  solution  of  (4.1)  on  f2x(-oo,oo), 
such  that  Uj^®0  on  Rx(-oo,0)  and  subject  to  the  boundary  conditions 

Uj^  =  uj  on  B^x(-cx),  00), 


on  B2x(-  00,  00), 


Here,  as  on  subsequent  occasions,  we  use  the  notation  for  the 
Stleltjes  Inverse  Introduced  In  Il](Theorem  1.3) •  Of .(2 .21), 

(2.22) . 
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which  follow  from  (2.9)  hy  virtue  of  (2.l),(2.3)f (2.6) .  On  the 
other  hand.  If  R  Is  simply  connected  and  (surface  tractions 

prescribed  over  the  entire  boundaxy)  the  unknown  stresses  are 
completely  characterized  by  (2.2),(4.3)*  and  the  second  of  the 
boundary  conditions  (2 .9) . 

Our  next  objective  Is  to  exhibit  a  particular  solution 
of  (4.2),  corresponding  to  vanishing  body  forces  and  to  a  given 
temperature  field  history,  by  generalizing  the  familiar  thermo- 
elastic  potential,  which  appears  to  have  been  discovered  original¬ 
ly  by  Borchardt  [12]. 

Ihenno-vlsooelastlc  potential .  Let  Op  (p*l,2)  and  0  meet  condi¬ 
tions  (a)  and  (b)  In  the  definition  of  a  viscoelastic  state. 
Suppose  0  l£  a  ( sufficiently  smooth)  real -valued  function  defined 
on  l2x(-oo,  oo)  that  vanishes  on  iex(-oo,o)  and  satisfies 

V^©  -  3a0*d{:^*d(2O^+G2)''^  (4.?) 

In  ftw(-oo,  oo) .  Then  throughout  that  domain  the  function  u  defined 

u  «  VO  (4.8) 

satisfies  (4.2),  provided  ^0  on  Rx(- oo,  oo) . 

To  confirm  this  assertion  one  need  merely  substitute 

from  (4.8)  Into  (4.2),  use  (4.7),  and  Invoke  the  algebra  and 

l8 

calculus  of  Stleltjes  convolutions.  A  suitably  regular  solution 
of  the  Poisson  equation  (4.7)  may.  In  turn,  be  generated  with  the 

T  y  T—  I  I  1  III 

For  related,  though  less  explicit,  particular  solutions  of 
(4.2)  see  Parkus  Il3](Chapter  Vl)  and  Nowackl  [7](Chapter  XI). 
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aid  of  a  Newtonian  potential,  as  Is  apparent  from  Lemma  9.1  in 
[1].  Indeed,  If  p  designates  the  right-hand  member  of  (4.7), 
such  a  solution  Is  given  by 


©(x,t) 


1  f  P(i#t) 

R  |x-L| 


for  all  (x,t)  In  Rx(-oo,  oo)  .(4.9) 


Ihe  last  theorem  evidently  supplies  yet  another  means  of  reducing 
the  boundary-value  problem  of  Section  2  to  a  mixed  problan  In  the 
Isothermal  theory.  Further,  the  present  theorem.  In  conjunction 
with  Theorem  9.2  In  [l],  leads  to  the  following  generadlzatlon  to 
thermo-vlscoelastlclty  theory  of  the  Papkovlch-Neuber  stress 
functions  In  classical  elastostatlcs . 

Generalized  Papkovlch-Neuber  solution .  Let  Gp  (p=l,2),0,©  meet 
the  hypotheses  of  the  preceding  theorem  and  let  P  conform  to  (b ) 
In  the  definition  of  a  viscoelastic  state .  Suppose  9  and  3!^^  are 
(sufficiently  amocth)  functions  defined  on  (?x(-oo,  00),  both  of 
which  vanish  on  <8x(-  oo,0)  and  satisfy 


with 


-  -  I  X’5  >  »  I  H  , 

H  «  P*dG^^*d(20^+02)"^  , 


(4.10) 

(4.11) 


In  f2x(-oo,cx3) .  Then  throughout  that  domain  the  function  u  defined 


u  B  70  +  v((p+x*t^)ed(0^+20g)  -  ^*d(20^+Q2)  (4.12) 


satisfies  (4.2) . 
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In  view  of  Theorem  9.^  In  [l],  the  foregoing  solution 
of  the  displacement  equations  of  equlllbrlxun  Is  complete  In  the 
sense  that  every  suitably  regular  solution  of  (4.2)  admits  the 
representation  (4. 7), (4.10), (4.12) .  Consequently,  the  Introduc¬ 
tion  of  the  displacement  potentials  (stress  fiinctlons)  ©,9,  emd 
reduces  the  mixed  bovindary- value  problem  to  the  determination  of 
appropriate  solutions  of  Poisson's  equation.  On  the  other  hand, 
since  9  and  ^  enter  (4.12)  imder  Stleltjes  convolutions  and 
because  of  the  structure  of  the  second  of  (4.6),  the  application 
of  the  bomdary  conditions  to  the  generating  stress  functions  will 
ordinarily  lead  to  a  system  of  simultaneous  Integral  equations. 

We  tiuti  now  to  the  well-known  correspondence  principle 
that  links  the  linear  theories  of  viscoelasticity  and  elasticity, 
nils  principle,  which  Is  of  paramovint  Importance  to  the  solution 
of  the  class  of  boundai'y-value  problems  with  which  we  ane  con¬ 
cerned,  may  —  with  the  aid  of  the  state  definitions  adopted  In 
Section  2  —  be  phrased  In  the  following  concise  form. 
Correspondence  principle .  Suppose 

~  Rx(-oo,  00)  (4.13) 

and 

^  B^x(-oo,  00),  =  S^  on  Bgx(-  00,  00) .  (4.l4) 

Assume  that  Qp(t)  (P®l,2),u^(x,t),e^j(x,t),P^(x,t),  and  0(x,t), 
for  every  x  In  R,  are  all  of  the  exponential  order  0(exp(s^t)) 
t  — >  00  ,  where  a  (real)  constant .  Then,  for  each  t)  with 

R«(n)  >  V 
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(4.15)^® 

where 

ii(n)  *  I  n^Cn)  »  >*(^)  “  3  n’SgCn)  *  (^.16) 

and  further 

Uj^(*jn)  “  “  ^(**n)  on  Bg.  (4.17) 

nie  theorem  just  stated  evidently  implies  that  if  the 
solution  to  the  thermo>vlscoelastlc  problem  governed  by  (4.13)/ 
(4.l4)  Is  Laplace-transformable/  It  must  coincide  with  the  Inverse 
transform  of  the  solution  to  the  one-parameter  family  of  statlonaiv 
thermoelastic  boxindary-value  problems  characterized  by  (4.13)/ 
(4.16),(4.17) .  Accordingly  the  correspondence  principle  yields 
a  reduction  of  the  original  problem  to  one  In  steady-state  thermo- 
elastlclty  theory. 

To  confirm  the  principle  one  need  merely  remove  the 
time -dependence  from  the  field  equations  and  boundary  conditions 
appropriate  to  a  viscoelastic  state  by  applying  the  Laplace  trans¬ 
form  to  (2.1), (2 .2)/ (2. 3)/ (2. 6)/ (2. 9)  and  compare  the  resulting 
systffln  of  aquations  with  (2,l)/(2.2),(2.3)/(2.l6),  and  (2.9).  In 
particular,  since  the  stress-strain  relations  (2.6),  because  of 
(2.8),  Imply  (2.12),  the  convolution  theorem  for  Laplace  transforms 
furnishes 

Recall  Footnote  No.  7  and  the  notation  for  the  Laplace  trans¬ 
form  Introduced  In  (2 .26) . 
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“  n'52(n)C^i5^( -^n)  -  3a§(*,Ti)], 


?  (».l8) 


Which,  by  virtue  of  (2.16),  account  for  the  "elastic  constants" 
(4.16). 


Suppose  now  that  the  formulation  of  the  original 
(thermo-vlscoelastlc)  boundary- value  problem  Is  based  on  the 
differential  operator  law  (2.30)  together  with  the  Inltiail  condl 
tlons  (2.32),  rather  than  on  the  relaxation  Integral  law  (2.6). 
In  this  case  the  conclusions  (4.15), (4.17)  continue  to  hold  true 
(in  the  presence  of  adeqmte  regiilarlty  assumptions  on  the  field 
histories  Involved),  provided  (4.16)  are  replaced  by 


as  Is  suggested  by  (2.33).  A  direct  proof  of  the  correspondsnce 
principle  applicable  to  dlfferentled  operator  laws  may  be  conducted 
by  means  of  an  argument  that  Is  strictly  analogous  to  the  one  we 
have  outlined  for  the  relaxation  Integral  law.  Note,  however, 
that  operating  on  (2.30)  with  the  Laplace  transform  one  arrives 
at 


Qg(n) 

w 


rej^(’,Tj)  -  3o0(*,fi)] 


>  (4.20) 


on  R,  for  all  t|  with  a  sufficiently  large  real  part,  only  If  the 
Initial  conditions  (2.32)  are  met.^®  Thus  (2.32),  the  physical 


For  details  see  the  px*oof  of  Theorem  4.7  In  [l]. 
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significance  of  which  was  discussed  In  Section  2,  are  precisely 
the  conditions  needed  to  Justify  the  usual  formal  application  of 
the  Laplace  transform  to  differential  operator  laws.  In  which 
Jump  discontinuities  In  their  relevant  time 

derivatives)  at  t®0  are  either  tacitly  ruled  out  or  else  (2,32) 
are  tacitly  asstimed  to  hold.  This  Issue  was  apparently  first 
recognized  by  Comellussen  and  Lee  [l4]  and  was  further  clarified 
by  Boley  and  Weiner  [3]  (Article  15.6);  It  Is  treated  rigorously 
In  considerable  detail  In  [l]. 

Tlie  preceding  analogy  between  boundary-value  problems 
In  the  linear  theories  of  viscoelasticity  and  elasticity  has  Its 
origins  In  a  paper  by  Alfrey  [15].  Alfrey's  form  of  the  corre¬ 
spondence  principle,  which  is  based  on  the  differential  operator 
law,  is  confined  to  Incompressible  solids  and  to  Isothermal  condi¬ 
tions.  A  limited  extension  of  Alfrey *b  analogy  to  compressible 
solids  Is  due  to  Tslen  [I6],  who  ab€uidoned  the  assumption  of 
Incompressibility  and  postulated  Instead  an  artificial 
(physically  unrealistic)  connection  between  the  devlatorlc  and 
the  dllatatlonal  response  of  the  material.  Integral  transforms 
were  first  brought  to  bear  on  the  question  raised  In  [I5]  by 
Read  [I?]*  who  reduced  the  viscoelastic  problem  for  a  compressible 
solid  to  an  elastostatlc  problem  with  the  aid  of  the  Fourier 
transform.  Althovigh  the  stress -strain  law  employed  in  [17] 
Involves  only  three  Independent  differential  operators,  this 
superfluous  restriction  does  not  affect  the  generality  of  Read's 
argument.  Brull  [I8]  arrived  at  an  analogous  correspondence 
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principle  for  a  medlvun  v/lth  a  continuous  spectrum  of  relaxation 
times  by  means  of  the  Laplace  transform  and  on  the  basis  of  the 
relaxation  integral  law;  in  [17 ]  the  volumetric  response  is 
(unnecessarily)  assumed  to  be  purely  elastic.  Lee  [19]  deduced 
the  comterpart  of  Brull's  version  of  the  correspondence  principle 
for  the  general  differential  operator  law.  Hilton,  Hassan,  and 
Russell  [20]  allowed  for  thermal  expansion  within  the  framework  of 
Alfrey‘s  analogy,  whereas  Lee‘s  scheme  was  similarly  extended 
to  thenno-vlscoelastlclty  theory  in  [21].  Additional  related 
restilts  will  be  referred  to  later  on.  We  note,  in  passing,  that 
the  initial  conditions  (2.32),  which  must  supplement  the  differ¬ 
ential  operator  law  (2.30),  are  not  mentioned  in  [15], [l6], [I7], 

[19 ]<  and  [20],  while  the  particular  initial  conditions  presupposed 
in  [21]  are  und\ily  stringent  and  hence  misleading. 

The  extended  coirespondence  principle  was  applied  in 
[21]  to  the  solution  of  a  particular  spatial  boundary-value 
problem  and  has  since  been  further  exploited  in  the  quasi-static 
analysis  of  thennal  stresses  and  deformations  for  linear  visco¬ 
elastic  solids  with  temperature-independent  properties .  Specific 
examples  and  pertinent  references  may  be  found  in  [3]/[7]»  and 
[13].  Nowackl's  [7]  monograph  supplies  a  convenient  guide  to 
recent  Polish  papers  on  problems  of  the  type  formulated  in 
Section  2.  Among  these  we  cite  investigations  by  Nowacki  [22], 

[23]  and  by  Sokolowskl  [24].  A  study  by  Shaffer  and  Levitsky  [25] 
should  also  be  referred  to  in  the  present  context.  Additional 

6t  ....... 

nie  inclusion  of  the  thermal  expansion  term  is,  strictly 
speaking.  Inconsistent  with  the  acc<»ipanylng  assumption  of 
IncOTipresslblllty . 
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studies,  which  aim  primarily  at  dynamic  thermal  effects  In 
temperature-independent  linear  viscoelastic  solids,  will  be 
mentioned  In  Section  7* 

5 .  Temperature-dependent  response ;  formulation  of  bovindary- 
value  problems  for  thermo-rheologlcally  simple  materials . 

It  Is  evident  from  the  preceding  section,  that  the 
quasi-static  analysis  of  thermal  stresses  and  deformations  within 
the  linear  theory  of  homogeneous.  Isotropic  viscoelastic  solids, 
and  In  the  absence  of  thermo-mechanical  coupling,  presents  no 
essential  difficulties,  provided  the  mechanical  response  of  the 
material  Is  assvimed  to  be  Independent  of  the  temperature.  Un- 
fortvinately,  as  was  pointed  out  In  the  Introduction,  such  a 
treatment  of  the  problem  Is  remote  from  physical  reality  unless 
the  temperature  range  Involved  Is  exceedingly  small.  The  remainder 
of  this  paper  Is  chiefly  concerned  with  the  modifications  arising 
In  the  theory  discussed  so  far  if  the  releucatlon  moduli  In  (2.6), 
the  creep  compliances  In  (2.20),  and  the  material  response 
parameters  In  (2.30),(2.3l)  are  themselves  temperatiu?e -dependent . 

The  Importance  of  effects  that  stem  from  the  temperature- 
dependence  of  the  viscoelastic  properties  was  emphasized  and 
Illustrated  relatively  early  by  Preudenthal  [26], [27], [28], [29] . 
Special  problems  concerning  temperature-dependent  linear  viscoelas¬ 
tic  solids  were  also  considered  by  Hilton,  Hassan,  and  Russell  [20], 
as  well  as  by  Hilton  [30]^  Rongved  [3l]#  Weiner  and  Mechanic  [32], 
Landau,  Weiner,  and  Zwlcky  [33],  and  Aggarwala  [34].  All  of  the 
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problems  treated  in  the  publications  Just  listed  concern  either 
infinite  plates  or  (solid  or  hollow)  spheres  and  clrcvilar  cylin¬ 
ders,  the  field  data  and  boundary  conditions  being  such  as  to 
preclude  the  presence  of  more  than  a  single  space  coordinate. 
Further,  the  \3nderlylng  stress-strain  relations  correspond  to 
degenerate  cases  of  the  differential  operator  law.  In  fact, 
the  solids  considered  exhibit  Maxwell,  Kelvin,  or  Standard 
Linear  behavior  in  shear,  while  the  volumetric  response  is 
ordinarily  taken  to  be  either  elastic  or  incompressible.  Finally, 
the  temperature-dependence  of  the  material  is  in  all  cases,  with 
the  exception  of  [3^]^  confined  to  the  viscosity  parameters  and 
its  specific  choice  is  usually  in  part  a  matter  of  analytical 
expediency.  Of  particular  interest  is  Rongved’s  [31]  explicit 
solution  (in  closed  elementary  form)  for  the  transient  thermeJL 
stresses  in  an  elastically  compressible  Meucwell  sphere  since  it 
accommodates  an  arbitrary  radially  symmetric  temperature  field 
history  and  presupposes  no  restrictions  upon  the  temperature-de¬ 
pendence  of  the  shear  viscosity. 

A  systematic  scheme  for  including  the  influence  of 
temperature  on  the  response  characteristics  of  a  viscoelastic 
solid  in  the  analysis  of  thermal  stresses  was  initiated  by 
Norland  and  Lee  [35  ]#  who  took  as  their  point  of  departure  the 
temperature-time  equivalence  hypothesis  originally  proposed  by 
Leaderman  [36]  and  subsequently  introduced  in  a  slightly  differ¬ 
ent  form  by  Ferrl  [37].  According  to  this  postvilate  the 


562(25)/19 


36 


mechanical  response  of  the  material  Is  affected  by  a  unlfom 
temperature  change  only  within  a  uniform  change  of  the  time  scale; 
the  response  Is  uniformly  accelerated  or  slowed  down  depending  on 
whether  the  tanperature  Is  Increased  or  lowered.  Materials  obey¬ 
ing  the  temperature -time  equivalence  postulate  are  said  to  be 
"thermo-rheologlcally  simple'*  In  the  terminology  of  Schwarzl  and 
Staverman  [38]  who.  In  their  expository  article  [39 cite  ea^erl- 
mental  evidence  In  support  of  the  postulate,  which  appears  to  be 
In  remarkably  good  agreement  with  tests  performed  on  a  variety  of 
high  polymers  over  considerable  temperature  ranges.  The 
analytical  framework  developed  In  [35]  was  further  explored  In 
[40]. 

We  proceed  now  to  a  brief  resume  of  the  theory  of  thermo- 
rheologlcally  simple  viscoelastic  solids  and  In  this  connection 
consider  first  the  required  generalization  of  the  relaxation 
Integred.  law  (2.6).  Thus,  let  GpCt)  (p=l,2)  henceforth  denote 
the  values  of  the  relaxation  moduli  at  time  t  measured  at  the  base 
temperattu?e  T^  and  denote  by^p(t,T)  the  corresponding  values 
measured  at  the  (fixed)  temperature  T.  Consequently 

Sp(t,To)  =  Op(t)  (p=l,2).  (5.1) 

The  temperature -time  equivalence  hypothesis  may  then  be  expressed 
analytically  by 

gp(t,T)  »  ap(5),  5  =  t9(T)  for  (t,T)  In  (-oo,oo)x[T^,T2],  (5.2) 

where  [T^,Tg]  Is  the  temperature  range  for  which  the  postulate  of 


562(25)A9 


37 


thermo -rheological ly  simple  behavior  Is  presumed  to  hold.  Here 
^  Is  the  "reduced  tlme"^  while  q>  represents  the  characteristic 
"shift  function"  of  the  material,  which  governs  the  contraction 
(or  the  stretching)  of  the  time  scale  that  results  from  the 
temperature  change  T-T^.  Clearly, 


9(T^)  -  1  ,  <p(T)  >0  (T^  <  T  <  Tg)  (5.3) 

and  9  Is  evidently  a  steadily  Increasing  function. 

The  constitutive  law  (2.6)  must  at  present  be  confined 
In  Its  validity  to  a  material  that  Is  permanently  maintained  at 
the  \mlform  temperature  T^.  In  these  circumstances  one  gathers 
from  (2.6)  on  setting  @«0,  recalling  (2.?)*  and  making  use  of  the 
commutativity  of  Stleltjes  convolutions,  that 


t  *=-oo 
t 


0^(t-t»)dej^j(x,t‘)i 


t  «=-oo 


^  (5.4) 


If  the  material  Is  Instead  permanently  at  the  vinlform  temperatiu*e 
T,  the  values  Op(t-t')  In  (5.4)  —  according  to  (5.2)  —  are  to 
be  replaced  by  Qp(5-^')>  where  lJ»«»t’9(T).  If,  finally,  the 
material  Is  iinder  the  Influence  of  a  variable  (time  and  position- 
dependent)  tenqjerature  distribution  T  with  values  In  [T^,T2], 

(5.4)  are  In  need  of  a  two-fold  additional  amendment:  first,  the 
definition  of  the  reduced  time  C  J»ust  now  be  generalized  to  account 
for  the  cumulative  effect  of  successive  temperature  changes; 
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second,  one  must  allow  once  again  for  the  Influence  of  thermed 

expansion.  Ihls  leads  to  the  modified  releucatlon  Integral  law 
t 


0ac(2.t) 


t  «»-oo 
t 

J  - 

t  '*-00 


t 

3a  f  02U-C«)de(x,t«), 
t‘i-00 


(5.5) 


with  ^ 

K  “  p(x,t)  «  J  <p(T(x,t‘))dt»,  =  p(x,t«).  (5.6)^^ 

0 


We  note  that  the  temperature  enters  (5.5)  both  through 
0  and  through  In  contrast  to  (2.6),  the  stress-strain 

relations  (5.5)  imply  a  nonlinear  dependence  of  the  local  instan¬ 
taneous  stresses  upon  the  local  temperature  history.  It  is 
apparent  from  (5. 6)^(5 .3)  that  p(x,  •)  is  a  monotone  increasing 
function  of  the  time  on  (-00,00),  the  Inverse  of  which  we  denote 
«(*,•) J  whence 


t  =  w(x,5) . 


(5.7) 


If  f  is  a  fxmctlon  of  position  and  time,  we  shall  consistently 
write  f  for  the  function  defined  by 

=  f(x,u)(x,5)) .  (5.8) 


- - 

Althovjgh  this  generalization  of  the  reduced  time,  introduced 
In  [35J^  is  convincing  on  physical  grounds,  a  rigorous  deduc¬ 
tion  of  (5.5)  from  (5.^)  and  the  temperature-time  equivalence 
postulate  appears  to  be  lacking. 
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By  subjecting  the  variable  of  integration  in  (5.5)  to  the  trans¬ 
formation  t*=w(x,5*)  one  may  eliminate  from  (5.5)  any  explicit 
dependence  upon  the  physlcail  time  in  favor  of  the  reduced  time. 
In  this  manner,  using  the  notations  adopted  in  (2.7)»(5.8)  and 
Involklng  once  more  the  commutativity  of  Stieltjes  convolutions, 
one  arrives  at  the  following  particularly  convenient  version  of 
the  modified  relaxation  Integral  law: 

“ij  *  L 

.  r  (5-9) 

■’kk  '  ('kk-3“®)*902  J 


As  is  apparent  from  (5.8)*  the  convolutions  in  (5.9)  are  taken 
with  respect  to  the  reduced  time  rather  than  the  physical  time. 

We  note  that  (5*9)  have  otherwise  the  same  structure  as  (2.6). 

Strictly  analogous  considerations  apply  to  the  generali¬ 
zation  of  (2 ,20)  and  ultimately  lead  to  the  modified  creep  Integral 
law  in  the  form 


y  (5.10) 


which  may  also  be  deduced  directly  from  (5.9).  Finally,  the 
same  process  of  specialization  that  leads  from  the  Integral  laws 
(2.6)  or  (2.20)  to  (2.30)  now  furnishes  the  modified  differential 
operator  law 

A 

Where  D  is  the  reduced-time  derivative  operator,  l.e. 


(5.11) 
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^  >  (5.12) 
while  the  polynomial  operators  Pp,Qp  (p«l,2)  retain  their  previous 
meaning.  To  (5.11)  one  must  adjoin  the  Initial  conditions  (2.32), 
which  remain  unaltered. 

On  referring  (5.11)  to  the  physical  time  t  by  means  of 
(5.6), (5.7), (5.8),  one  Is  evidently  led  to  a  pair  of  differential 
equations  which  has  the  same  structure  as  (2.30),  except  that  the 
previously  constant  response  parameters  PQ,„,qp  are  now  functions 
of  the  tfflnperature.  It  is  also  clear  that  these  functions  cannot 
be  prescribed  Independently  for  a  thermo-rheologlcally  simple 
viscoelastic  solid  since  the  entire  tonperature-dependence  of 
such  a  solid  Is  governed  by  the  single  shift  function  9.  More¬ 
over,  an  arbitrary  disposition  over  the  temperature-dependence  of 
the  response  parameters  in  (2. 30)  Is  not  only  Incompatible  with 
the  ten^erature-tlme  equivalence  hypothesis  but  is  also  readily 
seen  to  be  Inadmissible  on  energetic  grounds . 

The  Implications  of  the  equivalence  postxalate  for  elastic^ 
Maxwell,  and  Kelvin  solids  were  examined  In  [4o].  As  is  at  once 
apparent  from  (2 .15),  elastic  materials  with  temperature  depend¬ 
ent  moduli  do  not  belong  to  the  class  of  thermo-rheologlcally 
simple  viscoelastic  solids .  On  the  other  hand,  the  shear  response 
of  a  thenno-rheologlcally  simple  Maxwell  solid  Is  found  to  be 
characterized  by 

=ij  *  i  'ij  - 


(5.13) 
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In  which  n  and  t  are  the  shear  modulus  and  the  relaxation  time  of 
the  solid;  whereas  ^  remains  constant  in  the  present  instance^ 

T  obeys 

^  <  T  <  Tg)  ,  =  T(T^),  (5.14) 

and  is  thus  a  monotone  decreasing  function  of  the  temperature. 
Analogous  conclusions  apply  to  the  dilatatlonal  response  ajid  to 
solids  of  the  Kelvin  type.  It  is  interesting  to  note  that  the 
ad-hoc  assumptions  regarding  the  temperat\u?e -Independence  of 
certain  response  parameters  made  in  [26]  to  [33]  are  in  fact 
precise  consequences  of  the  temperature-time  equivalence 
hypothesis.  Before  leaving  the  present  subject  we  observe  that, 
as  pointed  out  by  Lee  and  Rogers  [4l],  the  shortcomings  of  the 
differential  operator  law  (i.e.  of  finite-spectrum  models)  are 
apt  to  be  even  more  pronounced  once  the  (thermo-rheologically 
simple)  temperature  dependence  of  the  material  is  taken  into 
account,  in  view  of  the  concomitant  contraction  of  the  time  scale. 

We  are  now  in  a  position  to  state  the  standard  mixed 
boxindary-value  problem  in  the  quasi-static  linear  theory  of 
thermo-rheologically  simple  viscoelastic  solids.  A  formulation 
based  on  the  modified  relaxation  Integral  law  may  be  phrased  as 
follows.  One  is  to  find  field  histories  which  —  for 

given  known  Q^,G2,a,T^,9,  and  prescribed  P^,T,u^,S^  — 

satisfy  (2.1), (2.2), (2. 3), (2.4), (5.5), (5. 6)  on  Rx(-oo ,oo )  and 
meet  the  initial  conditions  (2.8),  as  well  as  the  boundary  condi¬ 
tions  (2.9).  This  problem  suggests  the  subsequent  generalization 
of  the  definition  of  a  viscoelastic  state  Introduced  in  Section  2 . 
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Thermo-rheologloally  simple  viscoelastic  state.  We  say  that 
[ui,eij#Cij]  belongs  to  the  class  of  thermo-rheologlcally  simple 
viscoelastic  states  on  Rx(-oo,  oo)  corresponding  to  the  data 
0^,G2,a,T^,9,P^,T  for  the  temperature  rsmge  and  write 

[Ui,eij,Cij]^0l0j^,a2>«>T^#<P>Pi>T]  on  Rx(-oo,  oo)  (5.15) 


If: 


(a) 

(b) 

(c) 

(d) 


^1  g^aat  reqtAlrements  (a),(b)  In  the 
definition  of  a  viscoelastic  state ; 

<p  ^  continuous  on  ]  and  obeys  (5.3); 

T»T^  on  Rx(-oo,o),  ^  continuous  on  Rx[0,  oo),  and  has 
values  In  [T^jT^]; 

equations  (2.1), (2.2), (2.3), (2. 4), (5.5), (5. 6)  hold  on 
(2x(-oo,  00) , 


Evidently,  (5.15)  In  conjunction  with  q)(T)=l  for  T  In 
^^l'^2^  Implies  (2.10),  whence  this  specialization  of  the  shift 
fvinctlon  9  yields  a  reduction  of  the  present  theory  to  the  theory 
of  temperature-independent  materials  treated  In  the  preceding 
sections . 


6.  Implications  of  thermo-rheologlcally  simple  behavior, 
applications . 

Our  current  objective  Is  the  compilation  of  a  few 
general  conclusions  pertaining  to  the  theory  of  thermo-rheologlcal¬ 
ly  simple  viscoelastic  solids.  Most  of  the  theorems  about  to  be 
cited  are  extensions  of,  and  upon  the  appropriate  specialization 
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of  9  reduce  to,  analogous  propositions  In  Section  3.  We  state 
first  the 

Uniqueness  theorem.  Suppose 


[Ui,eij,c5ij]6  —  Rx(-oo,  oo) 

^  3^[0^,Q2,a,T^,(p,P^,T]  on  R^(-oo,  oo) 


r (6.1) 


and  let 


=  u^  on  B^x(  -  oo,  CO ) ,  =  S^  ^  B2^( "  ).  (6.2) 


Then 

[u^, e^ j,dj^ j ]  ®  [Wj^,0,0j  on  Rx(-oo,  oo),(6.3) 

where  w^=0  on  Rx(-oo,0)  and  represents  an  ( Infinitesimal )  rigid 
motion  of  the  entire  body  on  Rx[0,  oo) . 


This  result  Is  Implied  by  a  more  Inclusive  uniqueness 
theorem  established  In  [42],  the  latter  being  applicable  also  to 
ablating  viscoelastic  solids.  It  Is  Interesting  to  observe  that 
no  restrictions  on  the  shift  function  9,  beyond  continuity  euid 
(5.3)^  are  needed  to  assure  uniqueness  In  the  present  circumstances. 
The  next  theorem  may  be  Inferred  directly  from  the  definitions  of 
stationary  elastic  and  thermo-rheologlcally  simple  viscoelastic 
states  with  the  aid  of  (5.8),(5.9)*  and  Theorem  1.2  in  [l]. 

Initial  response.  Let  [Uj^,ej^j,Oj^j ]  meet  (5.15).  Then 
000  00 

[Ui,eij,<Jij]€£.[ji,)i,a,Pi,©]  on  R  (6.4) 


with 
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^  =  I  “  3  ®2' 

Consequently  the  initial  state  Is  again  a  stationary 

elastic  state  and  may  be  found  directly  from  the  initial  field 

and  surface  data.  We  turn  now  to  a  propostlon  which  follows  from 
2'^ 

a  result  cited  in  Section  3. 

Volume  change.  Let  ^  meet  (5*15)  and  suppose 

Qg  *  3Hh  on  ( -  oo,  oo  ),  (6.6) 

where  x  ^s  a  constant  (elastic  dllatatlonal  response) .  Then  the 
volume  change  is  given  by  (3.17)j(3.19) . 

Ihe  subsequent  restxLt  was  deduced  in  [6]  and  extends  an 
earlier  observation  due  to  Hilton  [9]. 

Stress -free  temperature  field  histories .  Assume  (5*15)  holds  and 
suppose  S^=0  an  Bxjo,  00),  P^s=0  on  Rx[0,  00) .  Then  on 

Rx(-oo,  00)  ^  and  only  if 

-  %(4)  +  ai^{t)x.^  for  (x,t)  Ijj  Rx[0, 00),  (6.7) 

where  a^^a^^  are  functions  continuous  on  [0,  00) . 

Related  theorems  concerning  stress-free  temperature 
fields,  applicable  to  plane  strain  and  generalized  plane  stress, 
are  given  in  [6],  where  the  two-dimensional  theory  of  thermo - 
rheologlcally  simple  viscoelastic  solids  is  explored  in  some 
detail.  The  two-dimensional  considerations  in  [6]  eilso  led  to  a 


See  the  remarks  following  (3.19). 
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connection  between  the  plane-strain  and  the  generalized  plane- 

oh. 

stress  solution  associated  with  one  and  the  same  plane  problem. 

The  chief  analytical  complications  Introduced  by  the 
departure  from  temperature-independent  viscoelastic  behavior  arise 
in  the  actual  solution  of  boundary-value  problems.  Uhforttihately, 
the  theory  of  integration  presented  in  Section  4  does  not  admit  a 
useful  extension  to  thermo-rheologlcally  simple  viscoelastic  media. 
This  is  true  in  particular  of  the  powerful  correspondence  principle 
discussed  at  the  end  of  Section  4,  which  yields  a  reduction  of  the 
viscoelasticity  problem  (for  a  solid  with  a  temperature -Independent 
response)  to  a  problem  in  steady-state  thermoelasticity  theory. 

To  bring  the  difficulty  to  which  we  have  just  alluded 
into  focus,  we  observe  first  that  the  integrals  appearing  in  the 
modified  relaxation  law  (5*5)  are  no  longer  of  the  convolution 
type.  Consequently,  an  application  of  the  Laplace  transform 
(with  respect  to  the  physlcetl  time)  to  (5.5)  falls  to  furnish 
algebraic  stress-strain  relations  in  the  transform  domain.  The 
alternative  version  (5.9)  of  the  modified  relaxation  law,  which 
has  the  requisite  convolution  structure,  suggests  the  possibility 
of  referring  also  the  remaining  field  equations,  as  well  as  the 
boundary  conditions,  to  the  new  Independent  variables  (x,^)  with 
a  view  toward  a  subsequent  elimination  of  the  reduced  time  by 
means  of  the  Laplace  transform.  Such  a  procedure  in  general  does 
not  result  in  a  worthwhile  simplification  of  the  problem.  Indeed, 

- - 

For  the  analogous  connection  in  two-dimensional  thermoelasticity 
theory,  see  Mindlln  sjid  Salvadorl  [43]#  p.762. 
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let  us  use  (5. 6), (5. 8)  to  effect  the  proposed  chemge  of  variables 
and  let  us  agree  to  write  f  .  and  f *  for  the  space  derivatives 
and  for  the  reduced-time  derivative  of  a  fvinctlon  f  of  position 
and  of  the  reduced  time.  Then  the  displacement-strain  relations 
(2.1)  become 


(6.8) 


whereas  the  stress  equations  of  equilibrium  (2.2)  pass  over  Into 


P 


=  0, 


(6.9) 


Because  of  the  terms  Involving  P  the  trainsforms  of  (6.8), (6.9) 

(taken  with  respect  to  the  reduced  time)  no  longer  possess  the 

desired  structure  of  (2.l),(2.2),  unless  p  .  vanishes.  Such  will 

be  the  case,  as  Is  clear  from  (5*6),  If  T  ^  a  function  of  the 

time  alone .  In  this  special  event  one  Is  accordingly  led  to  a 

significant  generalization  of  the  correspondence  principle.  To 

avoid  unduly  cumbersome  notation  In  the  subsequent  statement  of 

this  extension  we  shall  write 

00 

=  J  f(x,^  )exp(-Ti^)d5  (6.10) 

0 

for  the  Laplace  transform  with  respect  to  the  reduced  time  of  a 
function  f  with  veaues  f  (x,^  )  . 


Correspondence  principle  for  ptu?ely  time -dependent  temperature 
histories .  Suppose 
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with  T  a  function  of  the  time  alone,  and  let 

Ujj^  =  on  B^x(-oo,  oo),  ^  B2x(-oo,  oo)  .  (6.12) 

Assume  that  Op (C)  (p=l,2),  ^  T(^), 

for  every  x  R,  are  all  of  the  exponential  order  0(exp(s^5))  ^ 

^  >  oo ,  where  a  (real)  constant .  Then,  for  each  q  with 

ReCn)  >  s^, 

•,n)jeij(  •>Tl)i3^j(  •,T|)]€t  [ii(Ti),H(Ti),a,F^(  •,Tl),0(ri)  ]  on  R, 

(6.13) 

where 

ii(n)  =  §  ^\(n)>  >t(n)  “  ^  (6.i4) 

and  further 

Ui(  ’^n)  “  u5(  •,Ti)  on  B^,  S^(  .,n)  =  sj(  •,Ti)  on  Bg.  (6.15) 

The  preceding  analogy  becomes  trivial  If  B=B2,  Pj^=0  on 
Rx(-oo,  oo),  and  S^=0  on  Bx(-oo,oo).  In  this  Instance  the  analogy 
merely  confirms  our  previous  concluslon^^  that  Rx(-oo,  oo). 

In  view  of  the  well-known  theorem  on  stress-free  temperature  fields 
In  thermoelasticity  theory  ([3]j  Art.  3.9). 

Suppose  next  that  for  non-negative  time  the  temperature 
^  function  of  position  alone .  In  this  second  degenerate 
case  the  modified  relaxation  Integral  law  (5.5)  niay  once  again 
be  written  In  terras  of  convolutions  with  respect  to  the  physical 
time,  as  Is  apparent  from  (5.6).  In  fact  one  now  has 

Note  that  ©  conforms  to  (6.7)  since  T  Is  at  present  Independent 
of  position. 
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®lj  “  ^kk  “  (6.16) 

provided 

Hp(x,t)  =  Gp(t9(T(x)))  for  (x,t)  In  Rx[0, oo) .  (6.17) 

Upon  taking  the  Laplace  transform  of  (6.16)  one  finds  that 
Sij(x,n)  =  2^l(x,T^)e^J(x,n), 

'^lck(x^n)  =  3H(x,ti)rej^(x,Ti)-3aS(x,Ti)  J 

for  each  fixed  t)  with  a  s^lfflclently  large  real  part  and  for  all 
X  In  R,  if 

^•(x^Ti)  =  I  tiH^(x,n)>  x(x,T))  =-5n^^(x,n).  (6.19) 


^  (6.18) 


But  (6.18)  are  the  stress-strain  relations  appropriate  to  an 
Inhomogeneous  linear  elastic  solid.  Consequently,  if  the  tempera 
ture  depends  on  position  alone^  the  original  b oxindary- value 
roblem  In  thermo-vlscoelastlclty  theory  Is  reducible  to  a  stead; 
state  thermoelastic  problem  for  an  Inhomogeneous  medium.  Unhappily, 
this  generalization  of  the  correspondence  principle  Is  of  very 
limited  practical  Interest  since  the  complexity  of  the  reduced 
problem  Is  apt  to  be  comparable  to  that  of  the  original  one. 

For  the  differential  operator  law  an  analogous  extension 
of  the  correspondence  principle  to  temperature-dependent  materials 
was  given  earlier  by  Hilton  and  Russel  [44]  In  a  paper  which  Is 
based  on  their  prevloixs  report  [20]  with  Hassan.  The  analysis  in 
[44]  Is  confined  to  differential  operator  laws  with  temperature- 
dependent  material  parameters  and  thus  presupposes  a  medium  with 
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a  finite  relaxation  or  retardation  spectrum.  The  suialogy  deduced 
In  [44]  for  purely  time-dependent  temperature  distributions  Is 
limited  to  mechanically  Incompressible  solids  (although  thermal 
expansion  Is  included)  and  Is  an  extension  of  Alfrey*s  [15 ] 
principle.  On  the  other  hand,  no  restriction  Is  placed  In  [44] 
on  the  linear  viscoelastic  volumetric  response  for  the  case  of  a 
purely  position-dependent  temperature  field;  the  correspondence 
principle  obtained  In  this  instance  Is  an  elementary  generalization 
of  Read’s  [17 ]  work  In  the  isothermal  theory.  Finally,  [44]  con¬ 
tains  a  sketch  of  an  approximative  approach  to  the  general  case 
of  a  position  and  time-dependent  temperature  field  history,  based 
on  the  assumption  that  the  material  may  be  regarded  as  stepwise 
temperature-independent  (in  time) .  The  extent  to  which  such  an 
approximation  scheme  Is  computationally  feasible  remains  to  be 
assessed. 

We  turn  now  to  available  applications  of  the  quasl- 
statlc  linear  theory  of  thermo-rheologlcally  simple  viscoelastic 
solids.  As  was  mentioned  before,  all  but  one  of  the  speeded  In¬ 
vestigations  listed  at  the  beginning  of  Section  5  are  based  on 
temperature-dependent  viscoelastic  models  that  fall  within  this 
category,  eilthough  the  problems  treated  were  not  Intentionally 
selected  to  Illustrate  the  general  theory  under  present  consider¬ 
ation  . 


25 - 

In  connection  with  this  extension  of  Alfrey’s  scheme  the  authors 
remark  that  the  effect  of  body  forces  and  surface  tractions  may 
be  determined  separately  from,  and  may  be  superposed  upon,  the 
ensuing  thermal  stresses  and  deformations.  Such  a  superposition 
Is  In  fact  not  permissible  for  a  temperature-dependent  material. 
Indeed,  as  has  been  pointed  out  already,  a  purely  time-depend¬ 
ent  temperature  field  woxild.  In  the  absence  of  loads,  always 
Induce  vanishing  thermal  stresses. 
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TOie  temperature-time  equivalence  hypothesis  was  applied  - 
by  Norland  and  Lee  [35]  to  the  plane-strain  analysis  of  an 
Incompressible  hollow  circular  cylinder  which  Is  exposed  to  a 
stationary  radial  distribution  of  temperature  and  Is  subjected 
to  a  suddenly  applied  uniform  Internal  pressure. 

Implications  of  the  theory  In  circumstances  that  Involve 
temperature  distributions  depending  on  both  position  and  time  were 
studied  In  [40],  which  contains  exact  solutions  to  two  space 
problems  of  this  type.  Bie  first  of  these  concerns  an  infinite 
slab  which.  In  the  absence  of  loads,  is  confined  to  a  purely 
transverse  motion,  the  stress-inducing  temperature  field  history 
being  permitted  to  vary  arbitrarily  with  the  thickness  coordinate 
and  the  time.  The  second  problem  alms  at  the  thermal  stresses 
and  deformations  generated  In  a  sphere  by  an  arbitrary  transient 
radial  temperature  distribution.  Both  solutions  accommodate  an 
unrestricted  (thermo-rheologlcally  simple)  temperature  dependence 
of  the  material  and  apply  to  solids  with  a  continuous  relaxation 
spectrum,  except  that  In  the  second  problem  the  volimietrlc  response 
Is  assumed  to  be  elastic 

We  return  here  briefly  to  the  example  of  the  slab 
treated  In  [40]  since  certain  aspects  of  this  problem  have  a 
__ 

Although  this  assmptlon,  which  is  common  In  the  stress  analysis 
of  viscoelastic  solids.  Is  in  general  agreement  with  experimental 
findings,  the  available  quantitative  Information  regarding 
volumetric  viscosity  effects  appears  to  be  Inadequate. 
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wider  significance.  Thus  c^sider  an  infinite  plate  of  constant 
thickness  ,2a'  and  choose  the  coordinate  frame  in  such  a  way  that 
x^»0  cbtocides  with  the  middle  plane  of  the  plate.  Then  R  is 
the  region  characterized  by  -a  <  a  and  B  consists  of  the 
paif  of  bounding  planes  x-=ia.  We  now  seek  a  thermo-rheologicaaiy 
simple  viscoelastic  state^  ^  Rx(- 00,00)  correspond¬ 
ing  to  the  data  subject  to  the  bomdary  condi¬ 

tions 

^31  *  ®  Bx(-oo,  00).  (6.20) 

We  assume  further  that  T,  u^  are  functions  of  (x^^t)  alone  and 

«  Ui  ■  Ug  *»  0  on  Rx(-oo, 00).  (6.21) 

The  foregoing  constraining  assumptions  are  consistent  with  the 
governing  field  equations  and  boundary  conditions,  which  at 
present  are  readily  found  to  that  all  field  histories 

depend  solely  on  (x2,t),  while 

®11  “  ®22  “  ®ij  “  °  Rx(-oo,oo),  (6.22) 

^11  “  *^22'  **33  **1J  “  °  (I^J)  Rx(-oo,  00),  (6.23) 

whereas 

5U  ■' 

Note  that,  contrary  to  our  previous  asstmiptlon,  R  is  not 
bounded  in  this  Instance. 

See  [4oJ  for  details. 
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5  =  -e^KiQ^,  2^  =  (e-3a©)#d02  on  f2x (-00,00),  (6.24) 

If  we  conveniently  write  <5  and  e  for  and  respectively. 

The  elimination  of  e  among  the  pair  of  Integral  eqtia- 
tions  (6.24),  in  view  of  the  silgebra  of  Stleltjes  convolutions 
(Section  1  in  [l]),  now  leads  to 

^  A 

d  s=  -3aA*d©  on  J2x(-oo,oo  ),  (6.25) 


where  A  is  the  a\ixillary  response  function  defined  by 


A  =  Q^*dG2*d(20^+Q2)”^  on  (-00,00).  (6.26)^° 

Bearing  in  mind  (2.7), (5. 6), (5. 7), (5.8),  and  putting  x.^=x.f  we 

draw  from  (6.25)  that  the  desired  stress 

Integral  representation 
t 

d(x,t)  =  -3a  f  A(^-^  »)d©(x,t *)  for  (x,t)  in  Rx(-co,oo),  (6.27) 
t  00 


with  ^ 

K  =  p(x,t)  =  J  9(T(x,t'))dt‘,  =  p(x,t»).  (6.28) 

0 

Moreover,  (6.26),  by  virtue  of  Hieorems  1.2,  1.3  in  [1],  for  t 
In  [0,  00)  Is  equivalent  to 


o 

GA(t)  + 


0(t-t»)A(t«)dt«  »  L(t), 


(6.29) 


31 


where 


■55 - 

The  same  response  fvuaction  Is  also  found  to  play  a  crucial 

role  in  the  problem  of  the  sphere  [40]. 

Recall  (2.13). 
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G(t)  =  2a3^(t)+a2(t),  L(t) 


\\{t)  + 


0 

0 


(6.30) 


so  that  A  satisfies  a  linear  integral  equation  of  Volterra*s 
second  type.  Finally,  If  the  requisite  Laplace  trsinsfornis  exist. 


(6.29)  and  (6.30)  furnish 

^(n) 


2‘5^(Ti)+U2(n) 


(6.31) 


Formula  (6.27)  was  obtained  in  [40]  by  applying  the 
Laplace  transform  to  the  pair  of  integral  equations  (6,24)  and  by 
subsequently  inverting  the  solution  of  the  resulting  pair  of 
algebraic  equations.  As  was  noted  later  by  Lee  and  Rogers  [4l] 
and  Is  clear  from  the  preceding  derivation,  there  is  no  need  to 
Involve  the  transform  calculus.  A  similar  observation  is  relevant 
to  the  numerical  evaluation  of  the  solution  under  discussion, 
which  was  carried  out  In  [4o]  for  a  commercial  polyraethyl 
methacrylate  on  the  basis  of  available  relaxation  data  and  test 
results  for  the  shift  function  9.  In  this  connection  the  auxiliary 
response  function  A  was  determined  from  (6.31).  The  laborious 
procedure  adopted  in  [40]  consisted  in  first  devising  a  convenient 
analytical  approximation  to  the  relaxation  data  that  permitted  the 
analytical  computation  of  A;  the  required  valixes  of  A  were  then 
fomd  with  the  aid  of  two  asymptotic  inversions  —  one  applicable 
to  large,  the  other  to  small  values  of  the  time.  In  contrast,  Lee 
and  Rogers  [4l]  subsequently  re-compxited  A  by  solving  the  Integral 
equation  (6.29)  directly  on  an  electronic  computer.  Kielr 
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procedui^,  which  stays  closer  to  the  actually  measured  data,  is 
evidently  more  direct  and  also  more  accurate. 

The  values  of  A  obtained  by  the  two  alternative  methods 
described  above  are  compared  In  Figure  1,  In  which  E  Is  Yo\mg‘s 
modulus  of  the  Initial  elastic  response.  Figure  2,  which  Is 
taken  from  [4o],  displays  the  time -dependence  of  the  normal  stress 
d  In  the  middle  plane  of  the  slab  for  polymethyl  methacrylate. 

O 

In  the  \mderlylng  computations  it  was  assumed  that  a  /k  *  1  hour, 
k  being  the  thermal  dlffuslvity  of  the  material;  this  choice 
corresponds  to  a  slab  thickness  2a  of  approximately  5.7  cm.  The 
temperature  field  history  used  In  these  calculations  Is  that 
supplied  by  the  solution  to  an  elementary  transient  heat-conduc¬ 
tion  problem:  the  entire  slab  is  initially  at  the  uniform  base 
temperature  T^  «  80°C  and  Its  faces  at  t^O  are  suddenly  raised  to 
the  temperature  Tj^=110°C,  which  Is  steadily  maintained  thereafter. 

Curve  1  In  Figure  2  shows  the  stress  values  obtained 
from  (6.27) .  Curve  2  exhibits  the  behavior  predicted  If  the 
temperature-dependence  of  the  response  Is  disregarded  and  the 
analysis  Is  based  on  relaxation  data  appropriate  to  the  aversige 
surface  temperature  of  95°0,  Curve  3*  finally,  shows  the  results 
obtained  If  one  neglects  all  viscosity  effects  as  well  and  deter¬ 
mines  <J  on  the  supposition  of  purely  elastic  behavior,  using  the 
elastic  constants  characteristic  of  the  Initial  viscoelastic 
response.  As  Is  apparent,  the  three  curves  almost  coalesce  for 
an  initial  period  of  about  10  minutes,  during  which  the  "actual" 
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response  is  nearly  elastic  and  hence  temperature-independent. 

For  approximately  the  first  20  minutes  Curve  2  slightly  tinder- 
estimates  the  stress  values  predicted  by  Curve  1  since  during  this 
stage  the  temperature  of  the  middle  plane  is  less  than  the  average 
temperature  of  95°C  and  consequently  the  actual  relaxation  process 
is  slower  than  that  underlying  Curve  2 .  Later  on,  however. 

Curve  2  progressively  overestimates  the  stress  values  depicted  by 
Cuz^e  1  until,  after  a  two-hour  period,  the  actual  stress  is  less 
than  15  per  cent  of  the  corresponding  value  in  Curve  2.  These 
res\ilts  demonstrate  once  again  the  unrealistic  character  of  any 
viscoelastic  thermal-stress  analysis  that  leaves  out  of  account 
the  temperature-dependence  of  the  material  properties. 

The  relative  merits  of  the  direct  numerical  integration 
of  Integral  equations  over  the  use  of  Integral  transform  techniques 
In  viscoelastic  stress  analysis  were  further  discussed  and 
illvistrated  by  Lee  and  Rogers  In  [45].  Apart  from  the  fact  that 
this  approach  makes  more  efficient  use  of  experimental  data,  it 
obviates  the  need  for  long-time  extrapolations  from  the  data 
beyond  the  time  range  which  the  desired  solution  of  the  problem 
Is  to  cover.  Further,  such  a  treatment  extends  the  scope  of  the 
analysis  to  problems  which  are  not  amenable  to  a  Laplace-trans- 
fom  solution.  Among  these  are  mixed  problems  In  which  the 
surface  tractions  and  displacements  are  prescribed  over  time- 
dependent  subsets  of  the  boundary  and  problems  In  which  the 
boundary  Itself  Is  a  function  of  the  time  (siu?face  ablation). 
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An  example  of  the  latter  type  1b  dealt  with  numerically  In  [45], 
where  the  solution  to  the  problem  of  the  sphere  obtained  In  [40 ] 

Is  generalized  to  accommodate  an  ablating  spherical  body.  In 
addition,  [45]  contains  a  formal  solution  to  the  problem  of  an 
ablating  spherical  shell;  for  fixed  boundaries  this  problem  was 
also  treated  previously  In  [6], 

Although  the  developments  Just  described  are  promising. 
It  Is  well  to  keep  In  mind  that  boundary-value  problems  within 
the  theory  under  consideration  are  reducible  to  the  solution  of 
Independent  Integral  equations  only  In  highly  exceptional 
clrcximstances .  Furthermore,  all  of  the  applications  available  to 
date  depend  on  the  fortuitous  possibility  that  the  required  space 
and  time  Integrations  can  be  carried  out  separately.  While 
certain  less  degenerate  problems  may  nowadays  be  accessible  to 
a  ptirely  numerical,  treatment,  such  prospects  should  not  detract 
from  the  need  for  a  systematic  theory  of  Integration. 

7.  Concluding  remarks. 

Throughout  the  preceding  developments  the  material  has 
been  considered  to  be  Isotropic  with  regard  to  both  Its  mechanical 
and  Its  thermal  response.  For  an  anisotropic  linear  viscoelastic 
material  the  constitutive  relations  (2.6)  give  way  to 

‘’13  *  ('l3^13®>*^®13k.e' 

In  which  and  are  the  components  of  the  tensor  of 

relaxation  fxmctlons  and  of  the  thermal-expansion  tensor, 
respectively.  Furthermore,  and 
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nie  first  of  (7.2)  follows  from  the  symmetry  of  the  stress  tensor, 
while  the  second  entails  no  loss  In  generality  because  of  the 
symmetry  of  the  strain  tensor.  In  addition,  it  Is  usually 
assxmied  that 

“  ®kilj* 

Hils  last  set  of  symmetry  relations,  which  Is  Implied  by  (7.2)  in 
the  special  case  of  Isotropy,  expresses  an  Independent  reqxilrement 
Uie  available  theoretical  support  for  the  validity  of  (7.3) 
depends  on  thermodynamic  arguments  involving  an  appeal  to  Onsager* 
reciprocity  relations.^ 

If  (7.3)  are  presumed  to  hold,  the  generalization  to 
homogeneous  anisotropic  solids  of  most  of  the  theoretical  results 
discussed  in  this  paper,  presents  no  difficulties.  This  is  true, 
in  particular,  of  the  correspondence  principle  discussed  in 
Section  4  and  Section  6,  which  now  leads  to  a  tie  between  the 
linear  theories  of  anisotropic  viscoelastic  and  elastic  solids, 
as  suggested  by  Blot  [49].  An  extension  of  this  principle  to 
Inhomogeneous  (isotropic  or  anisotropic)  viscoelastic  solids  is 
equally  elementary, though  of  little  practical  consequence. 

op  '  '  ' ' 

^  See,  for  example.  Blot  [46], [47].  A  discussion  of  this  issue 
may  be  found  in  a  recent  note  by  Rogers  and  Pipkin  [48],  as 
well  as  in  [5].  Note  that  the  covmterpart  of  (7.3)  in 
elasticity  theory  is  implied  by  the  existence  of  an  elastic 
potential . 

In  this  connection  see  [20], [44],  as  well  as  Hilton  and 
Dong  [50]. 
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Next,  we  have  so  far  considered  only  the  quasl-statlc 
analysis  of  thermal  stresses  in  viscoelastic  solids.  (Hiermo- 
viscoelastic  inertia  effects  have  been  the  object  of  several 
recent  specisG.  investigations.  Among  these  we  cite  publications 
by  Katasonov  [51 3^  Nowacki  [23]»  [52],  [533#  [5^]i  and  Zorawski  [55 3^ 
[563, [573> [583.  The  physical  significance  of  these  investigations 
is  difficult  to  assess  since  they  invariably  presuppose  a  tempera¬ 
ture-independent  material,  do  not  Include  quantitative  discussions 

of  the  results  obtained,  and  usually  refer  only  to  Instantaneous 

-^4 

temperat\are  changes.'^ 

Further,  nothing  has  been  said  as  yet  about  thermo- 
mechanical  coupling  effects,  which  have  been  consistently  dis¬ 
regarded.  An  interesting  study  of  such  effects  is  due  to 
Hunter  [60],  whose  analysis  is  set  within  the  framework  of  the 
theory  of  thenno-rheologlcaLlly  simple  viscoelastic  solids.  As 
suggested  by  the  unassuming  title  of  [60 ],  this  subject  is  in 
need  of  fvirther  consideration. 

Finally,  it  may  be  well  to  emphasize  once  more  that  the 
present  paper  is  confined  essentially  to  linear  thermal  stress 
analysis  for  viscoelastic  solids.  Thus  no  attention  has  been 
given  to  the  Influence  of  nonlinear  viscosity  effects  and  to 
finite  deformations,  which  havint  the  realistic  treatment  of 
thermal  stresses  in  metals  at  elevated  temperatures. 

Acknowledgment .  The  author  is  greatly  Indebted  to  M.  E.  Ghirtln, 
who  read  the  manuscript  and  offered  numerous  helpful  criticisms 
and  suggestions. 

Related  studies  in  thermoelasticity  theory  (see,  for  example, 
[59]^  indicate  that  the  size  of  such  inertia  effects  is 
drastically  reduced  once  the  fiction  of  sudden  tanperature 
changes  is  abandoned. 
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Fig.  I  Aft )  for  polymethyl  methacrylate 
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©  Temperature-dependent  response:  80®C,  T,  *IIO®C 

©  Temperature -independent  response  based  on 
behavior  at  95®C 

©  Elastic  response  based  on  initial  behavior 


Fig.  2  Slob  problem.  Time-dependence  of  cr  in  middle 
plane  for  polymethyl  methocrylote 


